Confining Light in Deep Subwavelength Electromagnetic Cavities 
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We demonstrate that it is possible to confine electromagnetic radiation in cavities that are sig- 
nificantly smaller than the wavelength of the radiation it encapsulates. To this aim, we use the 
techniques of transformation optics. First, we present a "perfect cavity" of arbitrarily small size in 
which such confined modes can exist. Furthermore, we show that these eigenmodes have a contin- 
uous spectrum and that bending losses are absent, in contrast to what is observed in traditional 
microcavities. Finally, we introduce an alternative cavity configuration that is less sensitive to ma- 
terial imperfections and still exhibits deep subwavelength modes combined with high quality factor, 
even if considerable material losses are included. Such a cavity may be interesting for the storage 
of information in optical data processing and for applications in quantum optics. 



PACS numbers: 41.20.Jb, 42.70.-a, 42.79.-e 

Transformation optics has recently shed new light on 
the interaction between electromagnetic radiation and 
matter [IH3]. It provides a recipe to design components 
that guide electromagnetic waves along predetermined 
curved coordinate lines. The advantage of this technique 
is that it allows to approach an electromagnetic problem 
from a geometric perspective, by bending and squeezing 
the coordinate lines. These geometric distortions of space 
can then be converted into a medium with well-defined 
constitutive parameters Sj . Based on an early idea of 
Pendry [1] , transformation optics was first used to design 
a spherical perfect lens [5]. The most exciting example 
of transformation optics is the invisibility cloak, but it 
has also been applied for beam manipulation, lenses and 
illusion devices [2 [3] EHE] ■ The cloaking idea can more- 
over be used, e.g., in acoustics, to hide structures from 
acoustic waves [15 , hydrodynamics, to protect coastlines 
or platforms from tidal waves |16j , or quantum mechanics 
for cloaking of matter waves [17] . 

So far, the main focus of transformation optics has 
been on cloaking and beam manipulation. Here we want 
to show that these ideas can also be used to design de- 
vices that are able to confine electromagnetic energy. 
Nowadays, this can be achieved with microcavities — the 
most important implementations being Fabry-Perot, di- 
electric, and photonic crystal cavities [T51 [TOJ — and by 
the use of electromagnetically induced transparency to 
slow down or even stop light [20l HI] ■ The characteristics 
of a microcavity are determined by two important param- 
eters: the quality factor Q, which describes the temporal 
confinement of the electromagnetic field, and the mode 
volume V, which is a measure of its spatial extent [15] , 
Indeed, several applications involving optical storage re- 
quire electromagnetic energy to be confined in a small 
volume over a long period of time [18] . Unfortunately, 
traditional cavities are severely limited in size due to the 



wavelike nature of light, which imposes a lower limit on 
the mode volume and hence prevents the miniaturization 
of photonic components below the wavelength [55]. In 
addition, the electromagnetic storage systems mentioned 
above all suffer from fundamental losses, e.g., whispering 
gallery losses in dielectric microcavities [15] . 

In this letter, we want to present a dielectric cavity 
of deep subwavelength dimensions in combination with 
a high quality factor. We want to stress that we con- 
sider here a resonant cavity and not a waveguide struc- 
ture which always has free propagation — and hence no 
confinement — in one spatial direction. We start by re- 
calling the transformation-optical machinery leading to 
the invisibility cloak. The values of the permittivity and 
permeability that implement a distortion of the electro- 
magnetic space can be calculated by properly design- 
ing the transformation of the coordinate lines [3J [3]. 
In the case of an invisibility cloak, the electromagnetic 
fields cannot propagate inside the cloaked region, e.g., a 
sphere with radius R±. This is realized by mapping the 
physical radius R\ on the origin of the electromagnetic 
space. Additionally, the outer boundary at radius R2 in 
physical space is mapped onto itself in electromagnetic 
space, ensuring a smooth transition into the transforma- 
tion medium and eliminating reflections. Any continu- 
ous coordinate transformation r = /(r) that satisfies 
the boundary conditions f{R\) — and /(-R2) = R2 wm 
thus implement the effect of an invisibility cloak. 

For a cavity that encapsulates electromagnetic energy, 
we have to achieve the inverse of a cloak: light rays may 
not escape from the outer boundary for a cavity. We thus 
have to ensure that electromagnetic waves cannot pass 
beyond the outer radius i?2- Adapting the constraints 
used for the invisibility cloak, we can impose the bound- 
ary condition /(i?2) = at the outer boundary and al- 
low the radiation to penetrate into the inner region with 



2 




FIG. 1: The coordinate transformation of the perfect cav- 
ity. The surrounding space is made invisible through a ra- 
dial coordinate transformation that maps R2 on the origin 
in electromagnetic space and is matched with vacuum at Ri . 
Expressed in Cartesian coordinates (x l ), the coordinate lines 
get folded back on themselves in closed curves. This is the 
origin of the perfect confinement. 



radius R\. This amounts to imposing that the trans- 
formation function is continuous at this boundary. We 
therefore need the following boundary conditions for the 
transformation function of a cavity: 



f(R 1 ) = R 1 , f(R 2 )=0. 



(1) 



One can interpret this cavity as a medium that cloaks 
away the surrounding space, instead of the surrounded 
space. Let us first consider the cylindrical case as shown 
in Fig. 1. This setup consists of three regions: I and III 
are vacuum, whereas region II contains a transformation 
medium, where we will use a radial coordinate transfor- 
mation mapping the coordinates (p, <fr, z) onto the coor- 
dinates (p',(j>',z') as defined by the transformation 

Ri 



p 



{P-R2), 



(2) 



Ri — R-2 

while the other coordinates ((f), z) remain unchanged. 
This coordinate transformation, shown in Fig. 1, is one 
possible transformation that satisfies the boundary con- 
ditions derived above. As the radial coordinate gets 
folded, the corresponding coordinate lines in Cartesian 
coordinates follow closed loops: a particle, identified by 
a positive (negative) x-coordinate and moving along a 
vertical coordinate line in region (I), is bent towards the 
right (left) in the transformation medium and returns 
into the vacuum region on the same vertical line, staying 
bounded in this region for an infinite time. 

In order to confirm this geometrical picture, we have 
calculated the bounded modes of this system. Since the 
cavity is a linear system with cylindrical symmetry, the 
modes can be written as 



where m is an integer quantifying the angular momentum 
of the mode and w is an eigenfrequency. There is no z- 
dependence of the solutions, since we are considering here 
an infinite cylinder. In general, w can be a complex value 
u) = 10' + icj", where the imaginary part is related to the 
losses of the electromagnetic energy due to dissipation or 
radiation. The eigenfrequencies and their corresponding 
cigenmodes are determined by calculating the solutions of 
Maxwell's equations in each region (I, II, III) and combin- 
ing them with the proper boundary conditions. Without 
loss of generality we can assume TM polarization. In the 
vacuum regions (I) and (III), Maxwell's equations can be 
combined to Helmholtz' equation. One can derive that 
the magnetic field in these regions (I) and (III) is given 
by: 

Hf{p) = A J m (kop) + B Y m (k oP ), (4) 
Hf n (p) = E HgXkop) + F H$(k p), (5) 

where fco — w / c j Jm arL d Y m are the Bessel functions of 
the first and second kind, Hm^ and H$ are the Hankel 
functions of first and second kind, and (A,B,E,F) are 
complex integration constants; we may set B = F = to 
impose finite energy and Sommerfeld's radiation condi- 
tion. Inside the transformation medium (region II), our 
analysis has shown that the magnetic fields have a similar 
form with the radial coordinate replaced by f(p): 

H^p) = C J m (hf(p)) + D Y m (k f(p)), (6) 

where, once again, (C, D) are arbitrary complex num- 
bers. At the interface between two materials, the tan- 
gential components of the electric and magnetic fields 
must be continuous. When we apply these conditions at 
p = i?i and p = i?2, we find a set of four equations: 



A J m (k Ri) = C J m (k f{Ri)) + D Y m (kof(Ri)), (7) 
r f(Ri) 

f(Ri) v , 



A J'JkoRi) = C &lj> m ( ko f(Ri)) 

Hi 



D 



Ri 



-y^to/ORi)), 



(8) 



C J m {k f{R 2 )) + D Y m (kof(R 2 )) = E H$(k R 2 ), 



C ^JM(R2)) + D^Y^(k f(R 2 )) 
= EH'£\k Q R 2 ), 



(9) 



(10) 



E(r,t)=E(p)e- 



i(m<p—u)t) 



(3) 



where the prime (') denotes differentiation with respect 
to the radial coordinate p. 

Surprisingly, when we apply the transformation func- 
tion defined by Eqs. ([IJ, we notice that there is no quan- 
tization of the eigenfrequencies. This means that modes 
with an arbitrary value of ui' can exist inside this cav- 
ity, even if the free-space wavelength (Ao = 2ttc/uj') is 
many times larger than the dimensions of the cavity. In 
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(a) Cylindrical cavity: 
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FIG. 2: (a) Material parameters of the perfect cavity as de- 
fined by Eq. Q. (b) The parameters of the equivalent spher- 
ical implementation, (c) The magnetic field distribution of 
a perfectly confined cavity mode with R2/X0 = 0.0032; the 
magnetic field is exactly zero in the outer region, implying 
the absence of energy radiated away to infinity. 



(a) pIR, (b) xIR, 
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FIG. 3: The non-singular cavity is constructed from a perfect 
cavity, but with thin rims ARi and AR2 removed at the inner 
and outer boundaries (vertical broken lines), (a) The resulting 
material parameters do not assume extreme values, (b) The 
magnetic field distribution inside this cavity, corresponding to 
a deep subwavelength solution with R ou t/^o = 0.19. There 
is a small part of the mode situated outside the outer radius, 
corresponding to a quality factor Q = 1.1 X 10 10 . 



Fig. 2, we show the magnetic field of such a deep sub- 
wavelength mode in which the free space wavelength is 
three hundred times larger than the outer radius of the 
cavity. From this figure, one can observe that the field is 
exactly zero in the outside region (p > R 2 ). This means 
that the electromagnetic energy is entirely located inside 
the cavity: these subwavelength modes are thus charac- 
terized by an infinite quality factor. One might therefore 
call this device a "perfect cavity." Since the boundary 
conditions (|7|-(|To|) are independent of R\, the perfect 
cavity retains its properties for arbitrary values of the 
radius of the inner vacuum region. 

The values of the permittivity and permeability re- 
quired to materialize this cavity can be determined using 
the equivalence relations of transformation optics [21 [4] . 
With the transformation function as given in Eq. J2]), we 
find the following nontrivial components: 



P-R2 



p _ p _ P-R2 J, _ 6 _ P 

R\ p-R 2 
(fli - R 2 ) 2 P ' 



z\ = P z z = 



(11) 



(12) 



The variation of these components as a function of the 
physical coordinate p is shown in Fig. 2(a)-(b). Each 
component of e* and has a negative value, imposing 
the use of left-handed materials. It is easily seen that 
any transformation medium satisfying f{R\) — R\ and 
/(i? 2 ) = will have a region with left-handed materials. 
Subsequently, we notice the behavior at the outer bound- 
ary, which is analogous to the inner boundary of the invis- 
ibility cloaks: the radial component becomes zero, while 
the angular component tends to minus infinity. We also 
studied a spherical implementation of this cavity and we 
found that it also exhibits a continuum of deep subwave- 
length modes with perfect quality factor. The material 



implementation is more realistic than the cylindrical case 
as it does not require any component of the constitutive 
parameters that tend to infinity, as shown in Fig. 2(b). 

Unfortunately, we have found that the design discussed 
above is highly sensitive to the value of the material pa- 
rameters. Essentially being a strange kind of cloak, one 
might expect the same kind of sensitivity: as material 
parameters deviate from the ideal values, an invisibility 
cloak retains its cloaking characteristics, albeit less per- 
formant [BJ. In this case, however, when we perturb the 
cavity by taking away a little rim from the outer bound- 
ary, we notice that the eigenmodes disappear completely. 
In the second part of this letter, we therefore propose 
an alternative design of the cavity that eliminates this 
singularity. When the outer boundary is not perfectly 
mapped onto the origin in electromagnetic space, our 
simulations show that energy is radiated away to infin- 
ity. This prohibits the existence of confined modes. We 
can reintroduce deep subwavelength modes with an addi- 
tional perturbation at the inner boundary, thus removing 
the impedance matching. This will lead under certain 
conditions to destructive interference of the outside field. 
This idea is supported by the solutions of the dispersion 
relation and by our full-wave simulations with a finite- 
element solver (COMSOL Multiphysics) . We start from 
a perfect cavity with the transformation function 



Ri 



\JR\- R\ 



R\ 



(13) 



from which we now slice off thin rims at both the inner 
and outer boundaries, i.e., the material is situated be- 
tween the radii R\ + ARi and R2 — Ai?2 . The material 
parameters that constitute this transformation are shown 
in Fig. 3(a), where we have chosen perturbations of a few 
percent; this reduces the constraints on the materials sig- 
nificantly (the permittivity and permeability range from 
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FIG. 4: The magnetic field of an eigenmode of the cavity 
including material losses. When we fix the geometric param- 
eters at i?i = 0.00157? 2 , A-Ri = 0.757? 2 , AR 2 = 0.083i? 2 , 
a = 1.0, = 1.0 X 10~ 2 , and 7 = 1.5825 x 1(T 2 , we find a 
subwavelength solution with angular mode number m = 11 
at -Rout /A = 0.14 with Q = 1983. The real parts of e and fi 
stay bounded between and —5.2. 



—0.15 to —6.76). For this configuration, the dispersion 
relation allows for one single mode solution, for every in- 
teger value of the angular momentum parameter m. In 
Fig. 3(b), we plot the magnetic field of such a mode with 
to = 8, which has i?2/-V) = 0.19, i.e., the wavelength is 
more than five times larger than the outer radius of the 
cavity, and a quality factor of 1.1 x 10 10 . We notice that 
it has the same structure as for the perfect cavity. 

The reader might object that the typical high losses in 
the metamaterials required for the cavity's implementa- 
tion will adversely affect the subwavelength modes. We 
have therefore extended the formalism to calculate the 
eigenfrequencies of the cavity made from lossy metama- 
terials. By considering TM polarization, only one per- 
meability (/i 2 z ) and two permittivities (e p p , e ^ need to 
be implemented. We now introduce the loss tangents a 
and P of the permeability and the permittivity, respec- 
tively. We can show that the effect of these loss tan- 
gents boils down to replacing the transformation f(p) 
by g(p) = \JJ}- + icv)(l -I- ifi)f(p), and to add a factor 
(1 + i/3)- 1 to the right-hand side of Eqs. @-([lO]). We 
have determined the resulting quality factor by using re- 
alistic loss tangents that are found for fishnet structures 
operating at optical frequencies [23] . and that have been 
shown to converge well to bulk values for multiple-layer 
fishnets that we need here. The negative permeability, 
which occurs close to the magnetic resonance of the fish- 
net, has typically a large loss tangent of approximately 
a = 1, whereas the permittivities have a nonresonant 
Drude dispersion with typical loss tangent /3 = 10~ 2 . 
Our results show that such a lossy cavity still supports 
subwavelength modes, although with much lower qual- 
ity factor. We have found that the quality factor can be 
increased by improving the impedance matching at the 
inner boundary by filling the inner region with a nonmag- 
netic material with loss tangent 7 sa p. This procedure 



led us after optimization with COMSOL Multiphysics to 
the results shown in Fig. [4] we find again deep subwave- 
length eigenmodes (i? ou t/A = 0.14) with quality factors 
up to Q R3 2000. These Q-factors are higher than the the- 
oretical upper limit predicted for plasmonic cavities [21] ■ 

In this letter, we have demonstrated that electromag- 
netic radiation can be confined in cavities with dimen- 
sions well below the wavelength of the radiation. These 
cavities do not require materials with extremely high in- 
dices of refraction and therefore do not suffer from a 
breakdown of the effective medium approximation that 
would otherwise arise due to the discrete atomic nature 
of all real materials. 
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